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1 Introduction 


The past few years have seen much progress in our knowledge about the nucleon’s spin 
structure thanks to the experimental study of the spin asymmetries A^^x^Q'^) {N = 
p, n, d) in deep-inelastic scattering (DIS) with longitudinally polarized lepton beams and 
nucleon targets. Such experiments yield information on the ’longitudinally polarized’ (or 
’helicity-weighted’) parton distributions of the nucleon, which will be denoted as 

AlUx, q 2 ) ^ r{x, q 2 ) - r{x, Q^) , ( 1 ) 

where /^(x, (/^(x, Q^)) is the probability at scale Q of finding a parton-type / 

(/ = Qy Qy 9 ) ill ^ longitudinally polarized nucleon (for definiteness, say, a proton) with 
its spin aligned (anti-aligned) with the proton spin and carrying the fraction x of the 
proton’s momentum. Taking the sum instead of the difference in Eq. (^ one obtains the 
usual unpolarized parton distributions of the proton, 

/(x,Q^) =/^(x,Q2) +/^(x,Q2) . (2) 

Apart from /(x, Q^) and A 2 ,/(x, Q^), a third type of parton density can be defined 
which describes the properties of transversely polarized nucleons 11301.1: 

^rfiXyQ^) = - f^(x,Q^) , ( 3 ) 

where now /^(x, Q"^) (/'*'(x, Q^)) stands for the probability of finding a parton-type / in a 
transversely polarized proton with its spin parallel (antiparallel) to the proton spin. These 
parton densities are usually referred to as ’transversity’ distributions. They complete 
the twist-2 sector of nucleonic parton distributions and are therefore as interesting in 
principle as /(x, Q^) and Unlike the case of the latter two, there is no 

gluonic transversity distribution at leading twist p, 

Unfortunately, nothing is known as yet experimentally about the transversity dis¬ 
tributions. As was shown in i, S i, it is not possible to measure them in inclusive 
deep-inelastic scattering because they are chirally odd. It is widely accepted that Drell- 
Yan lepton pair production in collisions of transversely polarized protons offers the best 
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possibility to get access to the A^/ 0, 1^, ^ and the RHIC spin physics 

programme comprises experiments of this kind |TT[| . 


There being a shortage of experimental information on the transversity distribntions, 
there has been qnite some effort on the theoretical side, aiming at the qnestion of whether 
the Aif and A^f are entirely independent of each other. This is not expected really 
because one clearly has 

f^{x, Q^) + f^{x, Q^) = f{x, Q^) + f^{x, Q^) . (4) 


In fact, Soffer has derived an interesting non-trivial inequality [^, 


|At/|<-(/ + Az./) = /^, 


which is valid separately for each quark flavour. 


(5) 


Soffer’s inequality can be derived in the context of a parton model |^, 0. The 

question immediately rises of whether the inequality is maintained when QCD is applied 
T^, pA[|. There are two potential sources of concern here. Firstly, as indicated in 


Eqs. (0-(i. the parton distributions evolve with Q^. Since the corresponding evolution 
kernels at leading order (LO) [0,1 turn out to be different for the f{x, Q‘^), A^f^x, Q"^), 
ATf{x,Q^), it is not a priori clear that the inequality is preserved by QCD evolution. It 
turns out |^, however, that this is the case. The crucial point here is that the difference 
of the involved flavour non-singlet LO Altarelli-Parisi evolution kernels is always positive, 

i (p™(x) + AiP,™(x)) - ATq';>(x) = Cj.(l - x) > 0 , (6) 

where Cp = 4/3. In Eq. (|^) Pqq\ ApP^^ are the unpolarized and longitudinally polarized 
LO quark-to-quark splitting functions, respectively, which have been calculated in 


and are equal because of helicity conservation. ApPq^ is their transversely polarized 
counterpart of 0 . From Eq. (^) it follows that if one assumes validity of Soffer’s inequality 
(§) at some input scale Qq, the inequality will still be fulhlled at any higher scale, at least 
in the context of LO evolution. Strictly speaking, this statement only follows for flavour 
non-singlet, in particular valence, quark densities, which is due to the lack of gluons in 


the case of transversity distributions. However, as was shown in the picture does 


2 











not change when the full LO singlet evolution is taken into account for the f{x,Q‘^), 
ALf{x,Q^) on the right-hand side of (|^). 


Secondly, there is the question of what role higher order QCD corrections play in 
Soffer’s inequality. These enter in two ways, namely via the corrections to the splitting 
functions governing the evolution of the parton distributions, and as corrections to the 
cross sections for the processes which are used to extract information on the parton dis¬ 
tributions from experiment in order to see whether they satisfy Eq. As is well-known, 
the split-up of the higher order corrections into corrections to the splitting functions and 
to the cross sections is not unique but depends on the factorization scheme adopted, which 
in turn implies that the parton densities themselves are scheme-dependent beyond LO. 
Clearly, one is free to choose the factorization schemes independently for the f{x,Q‘^), 
^Lf{x,Q^) and Ar/(x, Q^), and it will in principle always be possible to violate Eq. (H) 
by an arbitrarily large amount by adopting certain ’sufficiently incompatible’ schemes for 
the three densities. One therefore has to ask in which factorization scheme(s) it makes 
sense to try to check Eq. beyond LO. In other words, what should be the meaning 
of Eq. in higher orders? This question has also been addressed in and will be 
one key issue of this paper. We will provide all ingredients for the next-to-leading or¬ 
der (NLO) framework of transversity distributions, which will serve to enable a check of 
Soffer’s inequality to NLO. 


As was pointed out earlier, the Drell-Yan process with transversely polarized hadrons 
appears to be the most likely candidate for measurements of the transversity distributions. 
Much theoretical work has therefore focused on this process pT[| . In par¬ 

ticular, the NLO corrections to the transversely polarized Drell-Yan have been calculated 
in and ^ within different factorization schemes. The result of the latter calcu¬ 
lation 0.0 was used in for a first assessment of the impact of NLO corrections on 
Soffer’s inequality. However, for future NLO analyses of (hopefully forthcoming) experi¬ 
mental data, knowledge of the NLO corrections to the Drell-Yan cross section alone cannot 
be sufficient as data will necessarily be taken at various different values of the invariant 
mass of the Drell-Yan muon pair which sets the hard scale for the process. According 
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to the previous paragraph, the NLO corrections to the Q^-evolution of the transversity 
distributions also have to be known in order to perform a complete and consistent NLO 
calculation. The main result of this paper will therefore be the presentation of the NLO 
evolution kernels for the transversity distributions - only when these are available can 
Soffer’s inequality be meaningfully studied beyond LO. The method we will use to cal¬ 
culate the NLO transversity splitting functions will be the one used by Curd, Furmanski 
and Petronzio in the unpolarized case |^, (see also @]) and which was also utilized 
recently to obtain the NLO evolution kernels for the longitudinally polarized parton 
distributions A^,/. 

The remainder of this paper is organized as follows. In section 2 we will outline the 
framework for our calculations. Also in our case the transversely polarized Drell-Yan 
process will be the starting point. Section 3 will give details of the calculation of the 
NLO transversity splitting functions. Section 4 will present the results and focus on the 
validity of Soffer’s inequality beyond LO. Finally, we will draw our conclusions in section 
5. 

2 Framework 

We begin by defining the spin-dependent cross section for dimuon production by protons 
with transverse polarization along, say, the a:-axis: 

rdATcr _ 1 frda^'' 
drdfj) 2 y drdcf) drdcf) J ’ 

where the superscript It (tJ.) denotes parallel (antiparallel) setting of the transverse spins 
of the incoming protons. Furthermore, r = M‘^/S, where M is the invariant mass of the 
muon pair and y/S the centre-of-mass energy of the hadronic collision; 0 is the azimuthal 
angle of one of the muons and is counted relative to the axis dehned by the transverse 
polarizations. For the discussion to follow in this section it is convenient to temporarily 
go to Mellin moment space where we can perform the Q^-evolutions more explicitly. The 
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Mellin moments of the cross section are defined by 


dcj) Jo drdcj) 


( 8 ) 


Including the NLO corrections to this cross section (as calculated in a certain factorization 
scheme such as the MS scheme), one has the generic expression^ |^, |] 
dAra^ ^ 


d(j) 


9S 


cos(20)Arid’^(M") 1 + 


27r 


-ArC?^’” 


where 




(9) 


( 10 ) 


Here we have restored for the moment the possibility of the scattering of two different 
hadrons labelled even though anything apart from pp scattering will presumably 

not be realistic experimentally. We will keep this notation when dealing with the function 
AxH, but otherwise suppress the indices in the following. The Mellin moments 

of the quark distributions Arq and the 0 {as) corrections AtC^^ are dehned in analogy 
with Eq. (I), e.g. Aj’g"'(Q^) = dxx^~^ATq(x, Q^). 

When including the 0{as) corrections in (P), it is crucial that the parton distributions 
are evolved to NLO accuracy as well. Since there are no gluons involved in the case of 
the transversity distributions, their evolution equations reduce to simple non-singlet type 
equations. Introducing 

Arql = Arq^ ± A^g" , (H) 

one has 


AtH'‘(Q^) = 


ATqlM)^Tqr2&) - ATqri{Q^)ATq’lM) 


( 12 ) 


and the evolution equations (see, e.g.. 


d 


dlnQ^ 



= ATP;,_(«.(Q^))ATg!!(Q^) , 

(13) 

{Arqf - Arq'Dm 

= A^P;,+(a,(Q2))(A^g!( - A^q'^m , 

(14) 



(15) 


^Here we only consider photon exchange. The contributions from exchange and interference 
can be straightforwardly included by making use of Eq. (3) of H. 
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where q, q' are different quark ffavours and Note that the first moment 

(n = 1) in Eq. (|I^) corresponds to the evolution of the nucleon’s tensor charge P, |^. The 
splitting functions /\TPqq,±{.(^s{Q‘^))i are taken to have the following perturbative 

expansion: 


= (?^) Arf + (^) At4‘’” + ■ ■ ■ . 


2ti 




2 ti 




(16) 


{a} = {gg, ±}, {EE}. As indicated, ^TPqq,^, ^rPqq- and AtPt.s are all equal at LO. 
It is convenient to write 


= AtP,<‘>’" ± AtP«>’" , (17) 

AtF^'” = ArPWf + ATF',',fe . (18) 

It will actually turn out that = 0, so that Eqs. ([r^),(p3D can be replaced by 

= ATP",+(a.(Q2))ArtI(Q2) . (19) 

Needless to say that the NLO splitting functions have to be known in the same factoriza¬ 
tion scheme as the corrections to the cross section, . 


The solution to Eqs. (0),(|T^ is well-known (see, e.g., ||23|) and reads 


= 1 + 


asjQl) - asiQ"^ 

ttPo 


2/3o 


X 




^TqliQl) 


( 20 ) 


MQDJ 

with the input distributions A 7 ’g^(Qg) at the input scale Qq. Here and in all previous 
NLO equations one has to use the two-loop expression for the strong coupling. 


a»(<J^ 


A lulnQVAIs 


4ir ftlnOVAfg ft’(•i„Q2/A|g)' 


( 21 ) 


with the QCD scale parameter A|^ and Pq = 11 — 2//3, Pi = 102 — 38//3, / being 
the number of active flavours. We note that all LO expressions are entailed in the above 
equations by setting the NLO quantities AtC^^^ AtP^}\ Pi, to zero. 


Formulas exactly analogous to the above hold for the unpolarized and the longitudi¬ 
nally polarized cases. The only complication here is that there are obviously contributions 
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by gluons to the Drell-Yan cross section and to the evolution of the quark singlet com¬ 
bination Let us assume for the moment that we are able to eliminate these 

experimentally by taking suitable non-singlet combinations of Drell-Yan cross sections for 
various conhgurations of scattering hadrons. In particular, when considering the differ¬ 
ence of the cross sections for pp and pp scattering, the function reduces to pure 

valence, 

elArci'l (Q") . (22) 

1 

where all parton distributions refer to the proton, and with similar expressions for the un¬ 
polarized and longitudinally polarized cases. We then only need to consider Eqs. (|9[), ([T3|) 
which exactly carry over to the unpolarized case (by omitting all At and replacing cos(20) 
by 2) or to the longitudinally polarized case (by replacing At by A^ and cos(20) by —2), 
respectively. 


We now adopt the Drell-Yan process as the process dehning the (NLO) parton distri¬ 
butions and follow |T^ to argue that (if at all) Soffer’s inequality for the valence densities 
should be valid beyond LO if we choose a factorization scheme in which the non-singlet 
Drell-Yan cross sections for the unpolarized, the longitudinally polarized and the trans¬ 
versely polarized cases all individually maintain their respective LO forms. This appears 
reasonable since the parton distributions are then directly related to a physical (observ¬ 
able) quantity. For example, in case of the cross section for transverse polarization it 
means 


dArcr'^ 




■cos(20)Ar^'^(M2 


(23) 


d(j) 9S 

where AtH^{Q‘^) is dehned as in ([T0|),(p^) but in terms of the NLO parton densities 
ATQ^iQ"^) in the new scheme. Inserting (pO]) into (^, expanding in and equating with 
(^) we hnd 


AtC(Q^)= 1 + 


«.(Qo) - 


TTpo 


AtS^ 


in the new scheme, where 


AtT!! = AtP^I^A - ^AtC. 


<?<?.- 


Po 

T 


' asm 

psm. 


-2ATPpf'"/l3o 


ATqfiQl) (24) 


DY,n 

1 


— AtP®’*" 

2/?o 


(25) 
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is a combination of (MS) NLO quantities which is manifestly independent of the factor¬ 
ization scheme due to the fact that the Drell-Yan cross section in Eq. (^) is a physical 
quantity. As a necessary condition for Eq. (|^) to hold for the valence densities, its mo¬ 
ments (for real n > 0) must also satisfy Soffer’s inequality. Thus 

2|ATr(Q')l 1 + («.(Qo) - a.(Q'))AT£:!^/(7r/?o) 


g!!(Q2) + Aigh(Q2) 


1 + - a^m) 




-2[/\T^qq -^qq 


29T(Qg) 


p(f )i^ 

TT rtrt 


X 


2\ATql{Ql)\ 


as{Ql 


q^m+^Lq^-m 


/(7r/3o) 
(26) 


has to be smaller than unity in the scheme we have dehned. Here we have used for the 
evolution of g!!(Q^)-|-A^g” (Q^) that and |^, ^ 

as well as |^, ^ where Pgq,±"^ = ^ Pqq’'^ is dehned in analogy 

with (1^). Clearly, the right-hand side in (^Bf ) will in general only be smaller than unity if 
this is the case for the input densities at Q = Qo, 2|ATg!!(Qo)l/(?-(i5o) + ALg!!(Qo)) A 1- 
This, of course, is impossible to prove or disprove within perturbative QCD. What one 
can do (and what is the purpose of this paper) is to assume validity of Soffer’s inequality 
for the input and see whether the inequality is preserved when going to higher Q^. From 
(I) we know that the exponent -2{ATP^f''^ - Pg®’”)//3o is positive, which demonstrates 
preservation of Soffer’s inequality at the LO level as discussed in [^. Even more, one 
hnds that for any reasonably small as{Qo) 


asiQl) - (^s{Q^) 

Tipo 


2f5o \ ^ 9" 




9fA p(^)’’ 

— 2[/AT-rqq 


' —P 

-L ( 


(0),r 


)//3o 


(27) 


remains 


smaller than unity. It also turns out that the unpolarized Pj^'^ is always negative. 


so that the corresponding term in the denominator of Eq. (^) will always help preserve 
Soffer’s inequality. What still is to be shown is that the remaining terms in (^) do 
not cause any problem if (in the ’worst’ case) 2|ATg!!(Qo)l/(^-(Qo) +‘^Lq-iQo)) = 1- 
Expanding in and taking into account the above observations one immediately hnds 
that in this case 




,n 


^rPqq,- 




DY,n 

q 




(28) 


has to be positive to guarantee validity of Soffer’s inequality for valence densities at NLO. 
Having worked out this condition in Mellin-n space for convenience, we will now return 






















to Bj0rken-a; space. In the following we will calculate the splitting function ATPql^-{x), 
to see whether the x-space counterpart of the expression in (^), 



is indeed positive. 

We emphasize that our considerations in the above equations only apply to the valence 
densities. We will also provide the splitting function which is needed for the NLO 

evolution of the combination ATq+ in (|T^. To examine, however, the validity of Soffer’s 
inequality for the non-valence quark densities would involve taking into account the singlet 
evolution of the unpolarized and longitudinally polarized parton distributions. This can 
only be done within a detailed numerical study, which is beyond the scope of this paper. 


3 Calculation of the NLO transversity splitting func¬ 
tions 


Our calculation of the NLO splitting functions for the transversity distributions closely 
follows the ones performed in the unpolarized ||T8|, p!9|, p0| and longitudinally polarized 


211 cases. We only briefly outline the method and its application to the transversity case 


here. A thorough overview of the technique can be found in |^, ^ 


The method we will use is set up in Bj0rken-x space. It is based on the factorization 
properties of mass singularities in the light-like axial gauge, specihed by introducing a 
light-like vector n {jp = 0) with n ■ A = 0, A being the gluon held. The general strategy 
then consists of hrst expanding the squared matrix element AtM for the scattering of 
two (transversely polarized) (anti)quarks into two ladders of two-particle irreducible (2PI) 
kernels [^. The crucial point is that in the light-cone gauge the 2PI kernels are hnite 
before the integrations over the sides of the ladders are performed. Collinear singularities 
therefore appear only when integrating over the lines connecting the rungs of the ladders 
This allows for systematically projecting out the singularities by introducing a pro¬ 
jector onto (transversely polarized) physical states, AtV. More precisely, AtV decouples 
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the product At{AB) of two successive 2PI kernels by projecting onto states of definite 
transverse polarization of the particle connecting the kernels and by setting this particle 
on-shell in kernel A. Transverse polarization for a quark entering a kernel is obtained by 
using 

u{p, s)u{p, s) = -^^75 , (30) 

for its spinor u{p,s). Here p is the quark’s momentum and s its transverse spin vector 
satisfying = —1, s-p = 0, s-n = 0, where the latter equality is a consequence of using 
n to dehne the longitudinal direction, 

n ■ p = pn ^ 0, n ■ t = p ■ t = 0 , (31) 

t being any transverse vector. The part of AtV that is relevant for the Dirac algebra is 
then found to be 

~ , (32) 

where k is the momentum of the particle emerging from the top of the kernel. As already 
mentioned several times, no gluonic projection operators are to be introduced in the 
transversity case. By means of AtV, AtM can be written in the factorized form 

AtM = At^ ® Aj'Ti ® Aj'r2 , (33) 

where Ayd is interpreted as the (hnite) short-distance (Drell-Yan) subprocess cross sec¬ 
tion, whereas the factors A^-T (one for each quark) contain all (and only) mass singularities 
and are process-independent. The explicit expression for AtT in terms of the 2PI kernels 
can be easily derived from |T^, Working in dimensional regularization (d = 4 — 2e 
dimensions) in the MS scheme it can be shown [^] that the residue of the 1/e pole of 
AtP corresponds to the evolution kernels we are looking for: 

ATr„ (j,a„i) (^) (1) , (34) 

where we have restored the subscripts {ii} = {qq, ±}, {HE} that distinguish between the 
various NLO quark-to-quark splitting functions. 

The Feynman graphs contributing to AtPm at NLO are obviously the same as in the 
unpolarized and the longitudinally polarized calculations |1l8| , |20| , |2T[| and need not be 
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repeated here. Let us instead collect a few important points of the calculation^: 


As mentioned above, the calculation relies on the use of the light-cone gauge. Fol¬ 
lowing 1^ ^ 1^ , we use the principal value prescription to regularize the spurious 
singularities resulting from the gauge propagator (see for details). 


From Eqs. (|30D,(|3^ we see that the Dirac matrix 75 enters the calculation via the 
projectors on physical states of transverse polarization of the quarks. In general, 
this can lead to complications when dimensional regularization is used. It turns out, 
however, that no problems related to 75 occur in our calculation due to the fact that 
all Dirac traces that we need contain two 75 matrices, one coming from the projector 
AtV in (|3^ , the other from (|30D . In this case, it is safe to use a fully anticommuting 


75 , which effectively removes the 75 from the traces via 7 I = 1. The only case where 
this does not work is the splitting function AtP^I'^ps for which there are two separate 
quark lines containing one 75 each. However, it is immediately obvious that for both 
Dirac traces the 75 is accompanied by an odd number of other Dirac matrices. So 
the traces vanish irrespective of the prescription for 75 . As a result, the splitting 
function ApP^^^ps vanishes identically. 


• The calculation is technically slightly more involved than in the unpolarized or the 
longitudinally polarized cases, owing to the appearance of scalar products of mo¬ 
menta with the transverse spin vector in the squared matrix elements. These scalar 
products, which are obviously always quadratic in s, introduce extra transverse 
degrees of freedom that need to be integrated out. In LO, this is a rather straight¬ 
forward task. The only scalar product of this kind that occurs is (fc • s)^, which after 
integration can be effectively replaced by 

(*: • 7 -- U ■ (35) 

This term is actually multiplied by e in the LO squared matrix element. As a 
result, it only contributes to the LO transversity splitting function in d = 4 — 2e 
dimensions, but not to its /o^xr-dimensional counterpart. In the 

^We use the program Tracer of for calculating the Dirac traces and performing contractions. 
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NLO calculation, however, the c?-dimensional LO result is needed. 

It turns out that the contribution from the term ~ e{k ■ s)^ is cancelled by other 
terms ~ e in the matrix element, so that for a; < 1 the d-dimensional LO splitting 
function is identical to the four-dimensional one of |^, 

0/y» 

ATP,';>-“‘(U = ArP™W=C'Fy^ (^^<1). (36) 

-L th/ 

Note that this result for the d-dimensional LO splitting function is at variance with 
the one of H in which residual terms ~ e{l — x) are present, presumably resulting 
from the assumption fc ■ s = 0 in |]^ . We will return to this point later. 


At NLO, the situation is more complicated. For the virtual graphs, besides terms 
involving [k ■ s) also powers of (r ■ s) appear, where r is the loop momentum. To 
integrate these, one needs tensorial two- and three-point functions. It turns out 
that all the functions we need here already appeared in our unpolarized 
longitudinally polarized calculations, so that there is no really new integral to 
be calculated. When calculating the squared matrix elements for the real diagrams, 
one encounters the terms {k ■ s)^, {k ■ s){li ■ s) and {li ■ s)^, where li is the momentum 
of one of the outgoing ’unobserved’ particles. One can integrate these terms in 
two steps: hrst, one integrates over the terms (/i ■ s), ■ s)^. These integrations 

have lengthy expressions, given in detail in Appendix A. After this integration, the 
squared matrix element still depends on [k ■ s)^ for which one hnds, similarly to 


^^ (fc2-2x(p-fc)) . (37) 

After this substitution, the resulting expression for the matrix element squared 
contains just terms that are familiar from the unpolarized calculation. It can then 


be integrated using the techniques developed in |^, 


For a transversely polarized anfiquark with a spinor v{p, s) one has expressions 
identical to (^),(^^. Thus - unlike the longitudinally polarized case - there is 
no extra minus sign when calculating the function AxPql^ ■ 


• The endpoint contributions to the transversity splitting functions, i.e. the contribu¬ 
tions ~ (5(1 — a:), are necessarily the same as in the unpolarized and longitudinally 
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polarized cases since they are just provided by the residue of the pole of the full 
quark propagator [|^, |2y, |2^ . 

We are now in a position to present the hnal results. 


4 Results 


For completeness, we begin with the full one-loop result 


2 x 3 


where the -|—prescription is dehned in the usual way: 

/ dzf{z){g{z))^= f dz{f{z) - f{l))g{z) . 

JO JO 

To write down our hnal result for the we introduce 

2 x 




{1-xU ’ 


= 


'TT^ dz ^ ,1 — Z, 

— In -, 


l + cc 


2 ; 


2 ; 


1 TT^ 

= — 2 Li 2 (—x) — 2 Inx ln(l-|-x) H—In^x-, 

2 6 


where Li 2 (x) is the dilogarithm. We then have (cf. Eq. (pT])) in the MS scheme 

AtF"±(^) = Aj-qpi) ± ArP,'!'(i) . 


where^ 


AtF£)(x) = cl 


1 — X — + 21n(l — x)^ Inx SrP^l^ix) 

+ - Y + 6^(3)') 5(1 -x) 


+ -^CpNc 


■(1 - i) + (y + F Inx + In^x - y ) StPS\^) 


(38) 


(39) 


(40) 


(41) 


(42) 


^Needless to say that the -I—prescription is not needed if the function multiplying the factor 1/(1 — a;)^ 
is vanishing at x = 1. Neither is it to be taken into account in the function x). 
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^+ ^-6C(3)V(l-^) 


12 


9 


Ic^Ti 


Inx- SrP^qKx) 


va«-'i 


AtP^I^x) = Cf(Cf-^Nc^ 


- (1 - x) + 2S2{x)6TP!jg\-x) 


where Cf = 4/3, Nc = 3, Tf = /Tr = //2 and C{3) ~ 1.202057. 


(43) 

(44) 


For checking Soffer’s ineqnality via Eq. we still need the 0{as) corrections to the 
short-distance snbprocess cross section for the transversely polarized Drell-Yan, {x), 

in the MS scheme. The corresponding resnlt was hrst presented in Q] where, however, 
dimensional reduction rather than dimensional regnlarization was nsed. The translation 
of the resnlt to dimensional regnlarization was provided in |]^. As a check of the ex¬ 
pression in 0, we employ an earlier resnlt for the transversely polarized Drell-Yan cross 
section at 0 {as) which was obtained in ^ by assnming off-shell glnons to regnlarize the 
appearing poles. Stndying in detail the strnctnre of the collinear singnlarities for both 
dimensional and off-shell regularization, it is straightforward to transform the result from 
one regularization scheme to the other. Starting from the result in |^, we obtain in this 
way for the MS scheme: 


ArCf'^ix] 


= C, 


8 x 


dn(l 


X 


Axlnx Qxhr X 


X 


X 


X 


-I- 4(1 — x) 


+ “®) '^(1 - 


(45) 


This result indeed coincides with the one in |]^ for the choice Axd = —5(1 — x) in that 
paper. On the other hand, the expression for A-pd that is claimed in [|^] to provide the 
link to the MS scheme (in dimensional regularization) is Ayd = —5(1 — x) -1-2(1 — x). The 
reason for this difference lies in the discrepancy between our calculation and the one of 


for the (4 —2e)-dimensional LO splitting function (see the remark after Eq. (^)). We note 
that the correctness of our result (^61) for this quantity is corroborated by the observation 
that the d-dimensional 2 —>■ 3 matrix element squared for the process qq (with 

transversely polarized incoming (anti)quarks) nicely factorizes into the product of the 
d-dimensional 2 —>■ 2 matrix element squared for gg —/i■*■/!“ times the splitting function 
of Eq. (|36D , when the collinear limit of the gluon to be parallel to one of the 
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incoming quarks is taken properly. This clearly demonstrates again the correctness of i 
The result of for the transversely polarized Drell-Yan cross section at NLO (obtained 
within dimensional regularization) therefore necessarily corresponds to a different (non- 
MS) factorization scheme. 


The remaining NLO ingredients for Eq. ( PU[ ) are the unpolarized quantities 
and which can be found in [^, ^ and |^], respectively. Inserting everything 

into Eq. (^) it turns out that all distributions at a: = 1 (-1—prescriptions and ^-functions) 
drop out, so that everything is regular at x = 1. We can therefore plot the result for 


D(x) = (q«(x) - - T ('^r(^) - ArC°''(i)) 


(46) 


directly in x-space. This is done in Fig. 1 for / = 3 active flavours. One can clearly see 
that D{x) is always positive, which demonstrates the validity of Soffer’s inequality for 
valence densities at the NLO level in the sense that if Soffer’s inequality is valid at the 
input scale, it will not be broken at any higher Q^. As becomes visible from Fig. 1, D{x) 
diverges at x —0. This feature is a consequence of a rather peculiar small-x behaviour 
of the NLO transversity splitting functions and subprocess cross sections: the quantities 
in Eqs. (|43|) , (]4^ ) , (^ ) all have a constant limiting behaviour at x —0. All terms that rise 
logarithmically in x are dampened by a factor x from the ’LO’ splitting term 5TPq^{x). 
This behaviour is in contrast to the one for Pjj^{x)^ C^'^{x) in (^6[), hence the rise of 
D{x) at small x. More precisely, we hnd from 


<? 

for X 


0 


AtPA(^) 


Cf i2CF - Nc) + 2Clx In^ x , 
—2Cf {Cf — Nc) x In^ x , 


(47) 

(48) 


that is, the constant terms even cancel in A 7 ’Pqg +(x). We note that the logarithmic 
terms in (|47D , (^8D can be recovered from Eq. (3.15) of where a resummation of small- 
X double-logarithms for the transversity densities was performed^. The constant terms 
in (|47|), however, were not anticipated to appear in ordinary Altarelli-Parisi evolution in 
3T|. The very mild small-x behaviour of Eqs. (|47|) , (|48|) will lead to a strong suppression of 


"‘Here one has to correct for a missing factor 1/2 in front of the term fyioj) in Eq. (3.15) of 1^ that 
got lost in the derivation of this equation from |Q. I am thankful to A. Vogt for helpful communications 
on this point. 
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the transversity distributions with respect to the unpolarized and longitudinally polarized 
ones during evolution at small x. 



Figure 1: The functions D{x), D{x) (see Eqs. (^) , (|4^) ) for / 


3 active flavours. 


As discussed in Sec. 2 we have chosen a factorization scheme in which the NLO Drell- 
Yan cross sections for the unpolarized, the longitudinally polarized and the transversely 
polarized cases all retain their respective LO forms. We have done this to deal with quan¬ 
tities that are manifestly scheme-independent. From a practical point of view, this may 
not be very useful since the unpolarized and longitudinally polarized parton densities are 
usually dehned in the MS scheme. In particular, when building models for the transver¬ 
sity densities, one would like to know whether Soffer’s inequality is also preserved under 
NLO evolution in the MS scheme. If this is the case then it appears likely to make sense 
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to assume validity of Soffer’s inequality also for the MS NLO input parton densities and 
not only for the parton densities dehned in our scheme. For this reason the dashed line 
in Fig. 1 shows only the quantity 

D(x) = Pfp(x) - ArF,'‘l(x) , (49) 

which corresponds to performing the evolution in the MS scheme. As can be seen, Soffer’s 
inequality will also be maintained here. 


For convenience, we present the Mellin-n moments of the transversity NLO quantities 
^T-Pgg^±, in Appendix B. These are useful for a numerical evaluation of our results. 

Fig. 2 shows the moments of the splitting functions versus real n > 1. As can be seen, 
the moments are all negative. Furthermore, the numerical contribution of the splitting 
function A^Pg^^’^ is very small for n > 1. For comparison we also show the moments of 
the LO transversity splitting function. The expression for the hrst (n = 1) moment of 
AxPgq^- turns out to be 


= —Cl - —CpNc + —CpTf 

8 ^ 72 ^ 18 


99 :- 


(50) 


This quantity participates in the NLO evolution of the nucleon’s tensor charge which is 
given by the hrst moment of the sum of the valence transversity densities |P]. 


5 Conclusions 


We have presented a calculation of the next-to-leading order splitting functions ATPq|±(x) 
for the evolution of transversity parton densities. The method we have used here is 
the one of p8|, [l8| which was applied to the unpolarized and longitudinally polarized cases 
IF, W| and [^|, respectively. We have used our results to examine the preservation 


m 


of Soffer’s inequality for valence quark densities under NLO evolution. It turned out 
that (in a suitably dehned factorization scheme) the inequality is indeed maintained at 
any higher if it is satished at the input scale, which provides strong support for its 
correctness. In other words, 0{as) corrections do not seem to invalidate the inequality. 
We note again that extending our results to non-valence quark densities ahords inclusion 
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Figure 2: The Mellin-n moments of the NLO transversity splitting functions AtP^I^- 
(solid line) and ^ (dashed line) versus real u > 1. For comparison the dash-dotted 

curve shows the LO result A^Pg®’"’. 


of the full singlet evolution of the unpolarized and longitudinally polarized densities and 
thus has to be subject to a detailed numerical study. Combining our hnding for the 
valence part with previous experience from LO indicates that no surprises concerning 
Soffer’s inequality are expected in the singlet sector. 


Note added 


After completing this work we received the papers |33, 341 in which the two-loop anomalous 


dimensions for the transversity distributions were calculated using the Operator Product 
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Expansion. The obtained resnlts shonld correspond to the NLO splitting fnnctions we 
have calcnlated. Indeed, onr results in Mellin-n space in Eq. (B.l) are exactly identical 
to those in and (after revision of that paper). We note that our results were hrst 
presented on the ‘Ringberg Workshop on High Energy Polarization Phenomena’, Schlofi 
Ringberg, Germany, 25-28 February 1997. 
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Appendix A 


In this appendix we provide the expressions for the integrations over the terms {li ■ s) and 
(Zi ■ s)^. Dehning 

A = 2 |^-(1-p) + y - (p-, (A.l) 

we hnd: 


(/i • .) 
(/i ■ 


A 


2 — 2x{p ■ k) 

1 


(/c -s) , 


l-2e 


2(^1 ■ p)- \- 


A2 


pn 4(/c2 — 2 x{p ■ k))^ 


A2 {k-sf 


4 (A;2 — 2 x{p ■ k)Y 

as the effective substitutions in the squared matrix element. 


notation of the momenta see 21 


(A.2) 

\ (k^s^ 

— 2x{p ■ k) 

(A.3) 

For kinematics and the 
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Appendix B 


In this appendix we present the Mellin-n moments of the NLO quantities j_(x), 

AtC^'^{x) (in the MS scheme) which are useful for a numerical evaluation of our results: 


qq^r] 


= Cl 


+ -CpNc 


17 1-r] 


134 22 

, , - —^i(n) + —S2{n) 

12 n{n + l) 9 ^ ^ 3 ^ ^ 


n 


+AS,{n) 2S,{n) - + 8S{n) - ^'(-) 


,n. 


-CpTf 

3 ^ 


■\ + fs,(n) - 2S2(n) 


= Cp 


,71, 


- 8S{n) + ^'(-) 


— —— + 4S'i(n) + 12 (5'3(n) — C(3)) — 8 + -tt^ 
n[n + 1) 3 


(B.l) 

(B.2) 


where rj = A. The sums appearing here are defined by 


Skin) 

^ W'K 

III 


s'kil) 

7=1 J 


~S{n) 

f-lV' 

1=1 ^ 

(B.3) 


Their analytic continuations to arbitrary Mellin-n (which depend on 77 ) can be found in 

3^. 
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